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Abstract 

We have found weak long range antiferromagnetic order in the quasi-two- 
dimensional insulating oxide KCr3(OD)6(S04)2 which contains Cr^+ S =3/2 
ions on a kagome lattice. In a sample with ~76% occupancy of the chromium 
sites the ordered moment is 1.1(3)//^ per chromium ion which is only one third 
of the Neel value gusS = 3fiB- The magnetic unit cell equals the chemical 
unit cell, a situation which is favored by inter-plane interactions. Gapless 
quantum spin-fluctuations (A/ks < 0.25K) with a bandwidth of 60K >> 
Tjv = 1.6K are the dominant contribution to the spin correlation function, 

S{Q,uj) in the ordered phase. 

PACS numbers: 71.30.-hh, 72.15.-v, 75.30.Kz 
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I. INTRODUCTION 



More often than not, local moments in insulating magnets acquire a finite expectation 
value, < S > at low temperatures. The only known exceptions are certain quasi-one- 
dimensional antiferromagnets (AFM) in which coupling between one dimensional units can- 
not overcome the tendency of each unit to form a cooperative quantum mechanical singlet. 
It is of great interest then to identify and study higher dimensional "moment free" insu- 
lating spin systems. We have made some progress towards this objective by identifying a 
quasi-two-dimensional magnet in which | < S > | is only 0.4(1) x S, where S is the spin 
quantum number of the magnetic ion. 

The material in question is KCr3(OD)6(S04)2, a quasi-two-dimensional AFM in which 
S =3/2 Cr^+ ions occupy weakly interacting kagome lattices. Quasi-two-dimensional AFM's 
containing spins on kagome lattices are likely to have reduced values of | < 5 > | as a 
consequence of the weak connectivity and triangular motif of the lattice. For Heisenberg 
spins for example it is known that | < S > \/S ^ 1 for S ^ oo and T ^ whereas 
I < S > I = for S* = 1/2 [|]. The critical value of S below which a moment free magnet 
exists at T=0 is however unknown 0. 

Only a few AFM kagome systems have been studied so far. SrCr9pGai2-9pOi9 is a quasi- 
two-dimensional oxide containing antiferromagnetically interacting S = 3/2 Cr^+ ions in 
kagome bilayers ||To|-[l4|. The material exhibits predominantly dynamic spin correlations 
although a spin-glass-like transition at ksTg/J ^ 0.03 causes condensation of weak two- 
dimensional static correlations with an ordered moment corresponding to | < 5 > | = 
0.63(5)5* and a correlation length of order twice the spin-spin separation [|14|-|1^. There are 
also recent reports of short range magnetic order in D30Fe3(OH)6(S04)2 |l3- In contrast the 
Fe-jarosites AFe3(OH)6(X04)2, with A=K, Na, Tl and X=S or Cr which realize a S =5/2 
kagome AFM all have been found to develop long range order at low temperatures with 
ordered moments | < S > \ ~ 0.75* [lT8| , p!9[ . 



In this paper we report experiments characterizing static and dynamic spin correlations 
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in the S =3/2 jarosite kagome lattice KCr3(OD)6(S04)2. Our principal result is that long 
range order similar to that found in the 5 = 5/2 jarosite systems does develop but the 
spin polarization, | < S* > |, is only 0.4(1) times the Neel value, S. Quantum fluctuations 
account for the remaining spin of the Cr^"*" ions. The bandwidth of the fluctuation spectrum 
is approximately 6 meV and we place an upper bound of 0.03 meV on any gap in the 
spectrum. 

II. EXPERIMENTAL DETAILS 

The neutron scattering measurements were carried out on thermal neutron triple axis 
spectrometers at the NIST reactor. Measurements were taken with flxed incident neutron 
energies of 13.7meV or 5meV using PG or Be as fllters in the incident beam. Because impor- 
tant conclusions of this paper come from absolute measurements of the neutron scattering 
cross section, details of how we normalize scattering intensities from a powder sample are 
given in appendix A. 

A 2.8g powder sample of KCr3(OH)6(S04)2 and a 6g powder sample of KCr3(OD)6(S04)2 
were synthesized using previously published methods . We denote these samples A and B 
respectively. Powder neutron scattering conflrmed that the samples were single phase with 
the hexagonal alunite structure and low temperature lattice parameters a=7.2371 A and 
c=16.9544 A. Fig. |I] shows T = 20 K diffraction data from sample B with the Rietveld flt su- 
perimposed. In sample A the Cr^+ occupancy was 100% as determined from the slope of the 
inverse high temperature Curie- Weiss susceptibility versus T. In sample B prompt gamma 
activation analysis, reflnement of powder neutron scattering (see Table 1 and Fig. [l|) and 
chemical analysis indicated Cr^^ occupancy of 76(5)% whereas high temperature susceptibil- 
ity measurement indicated 90% Cr occupancy. Unfortunately problems with stoichiometry 
are not uncommon in this complicated class of materials All neutron scattering exper- 
iments reported here were performed on the deuterated sample B. 

The chromium lattice in KCr3(OD)6(S04)2 is illustrated in Fig A formula unit contains 
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three kagome planes of chromium ions stacked along the c-axis. Neighboring kagome planes 
(one shown with solid lines, another with dashed lines) are displaced by (|, |, |). A, B and 
C indicate three sublattices of the magnetic structure to be discussed later. Cr'^+ ions are 
coordinated by six oxygen atoms located at the vertices of trigonally distorted octahedra 
T8| , |22[] . Neighboring octahedra share one oxygen atom which is expected to mediate the 



strongest exchange interaction in the system between neighboring chromium ions separated 
hj d = 3.624 A in a kagome plane. Susceptibility measurements on samples A and B 
showed Curie-Weiss temperatures, Qcw, of -70(5)K and -54(2)K respectively. If nearest 
neighbor exchange in the kagome lattices indeed dominate, this corresponds to a bond 
strength J = {3kB\Qcw\/ ^S{S + 1)) = 1.2{l)meV where z = Aricr/Q is the average nearest 
neighbor coordination {ncr is the number of chromium atoms per unit cell determined from 
Rietveld analysis and listed in Table 1). Other potentially significant exchange interactions 
in order of the strengths expected are the interplane coupling (Cr-K-Cr, d=6.025(l)A) and 
the second nearest intraplane interaction (Cr-O-K-O-Cr, d=6.268(l)A). 



III. EXPERIMENTAL RESULTS 

A. Antiferromagnetic long-range-order 

Fig. ^ shows the temperature dependence of magnetic neutron diffraction at Q = (012) 
and specific heat, C/T, which reveal a magnetic phase transition at T/v = 1.8 and 1.55K 
<^ Qcw for samples A and B respectively. The entropy change between 0.2K and 2.9K is 
only 0.15-Rln(2S+l) and 0.16Mn(2S+l) respectively indicating that the ordered moment in 
these samples may be small and/or there may be substantial short range order above the Neel 
temperature. Note that in contrast to D30Fe3(OH)6(S04)2 [|17| where impurities induce long 
range magnetic order, it is our cleanest sample. A, which has the most pronounced specific 
heat anomaly at T/y. This indicates that long range order is the clean limit low T state of 
KCr3(OD)6(S04)2. 
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Fig. ^ shows low temperature magnetic diffraction data for sample B. The data was 
taken without final energy analysis and with an incident neutron energy, Ei = 13.7 meV, 
which far exceeds the energy range with appreciable inelastic magnetic scattering. The data 
therefore probes the equal time spin correlation function S{Q) = J^^S{Q,Lj)hduj. The two 
frames show low temperature (T=0.4K) data with two different backgrounds subtracted. 
Using diffraction data from immediately above the Neel temperature (T=3.6K) as a back- 
ground yields simply resolution limited magnetic Bragg diffraction indicating that the low 
temperature phase has AFM order. From comparing the line shapes of magnetic and nuclear 
Bragg peaks we conclude that the AFM correlation length exceeds 500 A. Magnetic Bragg 
peaks occur at the same values of Q as nuclear Bragg peaks which implies that the chemical 
unit cell described above is also a unit cell for the magnetic structure. This implies that 
all triangles of spins which are connected by in-plane Bravais vectors of the crystal have 
equivalent spin orientations which we denote by A, B, and C. To determine the actual spin 
directions and the relative orientation of spins in neighboring kagome planes we compared 
the measured Bragg intensities to those calculated for specific model structures. Table 2 
makes the comparison to the spin configuration shown in Fig. |^ (a). In this structure, spins 
labeled by A, B, and C were taken to lie within the basal plane and to be rotated 120° 
with respect to one another. Subsequent kagome layers are oriented so that each interplane 
triangle of chromium ions, such as those shaded in Fig. H, contains all three spin orientations 
: A, B, and C. This structure provides the best account of the measured intensities out of 72 
high symmetry spin configurations which were examined. The ordered moment determined 
from this analysis was \M\ = 1.1{3) fiB per occupied chromium site which if g ^ 2 and 
S = 3/2 implies that | < S* > | = 0.4(1)5'. Note that the error bar quoted here includes 
the uncertainty associated with other choices of the interlayer registry which also provide a 
reasonable account of the diffraction data. 
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B. Equal-time short-range spin correlations 



Since the long range ordered moment is small, sum-rules imply that we should find 
static or dynamic short range order which accounts for the remaining spin on Cr^"*" ions. This 
aspect of the problem is examined in Fig. ^ (b) which shows diffraction data similar to those 
discussed above, though this time using T=50KRi Qcw data for background subtraction. At 
T=50K we expect that magnetic scattering is almost Q-independent in the limited Q-range of 
the experiment. This is also the case for the thermally activated inelastic phonon scattering 
which is present in the high temperature data as well [^. In Fig. ^ (b) these terms are 
over-subtracted but the Q-dependence of the difference data should follow that of the low 
temperature equal time correlation function S{Q). For the purpose of comparisons to theory 
we determined the base-line for S{Q) by requiring that the data (including magnetic Bragg 
peaks) satisfy the total moment sum rule 3 x / dQQ^S{Q)/ J dQQ^ = ncr{g/'2YS{S + 1) 
using ncr = 6.9 as derived from powder diffraction. The most striking feature of the data is 
the unusually small fraction of the equal time correlation function which is associated with 
magnetic Bragg peaks. The result indicates that most of the chromium spin in this material 
remain spatially disordered down to T=0.4K. 

C. Energy resolved magnetic neutron scattering 

It is important to determine whether the short-range-order is static or dynamic. To 
settle this question we used the full triple axis configuration of the neutron spectrometer to 
measure the low temperature (T=70mK) spectrum of magnetic fiuctuations at two values 
of Q. The result is shown in Fig. |^. Contributions to the detector count rate not associated 
with inelastic scattering from the sample were measured and subtracted as described in the 
figure caption. Integrating the data over hu we obtained the values for S{Q) shown with 
solid symbols in Fig. ^ (b). Inelastic scattering with 0.4meV< hu < lOmeV clearly accounts 
for most of the equal time correlation function and hence we may conclude that the short 
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range order in KCr3(OD)6 (804)2 is mostly dynamic in origin. From the energy spectra 
shown in Fig. ^ we estimate a magnetic bandwidth of approximately 6meV. In addition we 
can place an upper limit of 0.4meV on any gap, A, in the magnetic excitation spectrum. A 
stricter constraint on A is derived from the specific heat data of Fig. ^ where the absence 
of activated behavior for T > 0.25 K indicates that A > 0.03 meV. 

Further information on magnetic fluctuations is provided by the constant-Zici; scans shown 
in Fig. 1^. The data display an asymmetric peak centered close to the strongest AFM 
Bragg point, r = (101). The peak progressively broadens as energy transfer is increased. 
The sharp leading edge and gradual decrease in intensity with increasing Q is indicative 



of low dimensional magnetic correlations [|^. Based on inelastic constant-Q and constant- 
huj scans we can also evaluate the a lower hound on the total fluctuating moment 3 x 
/ S{Q,huj)Q'^dQhduj/ J Q^dQ = 15 per unit cell. This number is an appreciable fraction of 
the total moment sum ncr{g/'^)'^S{S + 1) 26 per unit cell which indicates that the spin 
system remains predominantly dynamic at low temperatures in this material. 

IV. DISCUSSION 

A. Magnetic structure 

Having presented the experimental results we now turn to comparison with relevant 
theory. The type of magnetic order in KCr3(OD)6(S04)2 is surprising. Quantum as well 
as classical theories of the kagome magnet predict that the so-called x \/3 in-plane 
structure [Fig. ^ (b)] is favored over the "q=0" structure which we observe. In addition 
another kagome related material SrCr9pGai2~9pOi9 shows short range x -\/3 type order 



T^,[T6|. It is therefore surprising to find q=0 type AFM order in KCr3(OD)6(S04)2. We 
believe that interplane coupling plays an important role in stabilizing this type of order 
in jarosite systems. As is apparent from Fig. |^ (a), q=0 kagome AFM layers can be 
stacked to take advantage of ferro- as well as antiferromagnetic nearest neighbor interlayer 
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coupling, but -y/S x -^/S layers cannot. The reason is that the (|, I) in-plane translation 
between neighboring kagome planes is incompatible with the periodicity of the x 
magnetic structure and it is therefore impossible to arrange for all interplane triangles of 
spins to contain all three spin orientations A, B, and C as required to satisfy interlayer 
coupling. Indeed, all jarosite kagome systems studied so far have q=0 type in-plane AFM 
order. Interplane coupling also controls another aspect of the three dimensional magnetic 
structure : the length of the primitive magnetic unit cell in the c direction: Three identical 
but shifted q=0 layers can be stacked to satisfy AFM interplane coupling. For ferromagnetic 
interplane coupling, however, it is favorable to include time reversed copies of each of these 
three layers in the stack which leads to a six layer sequence and a doubling of the unit cell 
in the c-direction. Such unit cell doubling is actually observed for KFe3(OH)6(S04)2 ||T8[| . 



B. Comparison to spin wave theory 

We have compared our inelastic data to specific models for dynamic correlations in the 
kagome AFM. Several authors have considered spin waves in a long ranged ordered classical 
kagome magnet and there are a few unusual features to note. As a consequence of a local 
continuous degree of freedom (the weather vane mode) spin wave theories to lowest order 
in 1/5 yield a zero energy mode of excitations for all Q. Furthermore it can be shown that 
all classical long range ordered ground states lead to the same dispersion relation for the 
finite energy modes. We used the results of Harris et al. to calculate the inelastic neutron 



scattering cross section associated with the finite energy modes of the =0 structure p5[ . The 
calculated cross section was convoluted with the experimental energy resolution function 
and is shown as dashed lines in Fig. ^ The comparison involves only a single adjustable 
parameter which is the overall scale factor. The exchange constant was fixed to ImeV, the 
value extracted from high temperature susceptibility measurements. Simple lowest order 
two dimensional spin wave theory clearly cannot account well for the data. Specifically the 
sharp features predicted by the theory are absent in the data. This is actually quite unusual 
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for two dimensional magnets. The qualitative features of spin dynamics in La2Cu04 are for 
example very well described by such theories [ p^ . 

C. Equal-time spin correlations 

Theoretical predictions are also available for the equal time correlation function S{Q). 
We compare these to our data in Fig. ^. The dotted line is the result of a Monte Carlo 
study |2^, and the dashed line resulted from a high-T expansion p|. Given that the ex- 
perimental data probe low temperature spin correlations it is surprising that apart from the 
magnetic Bragg peaks the high temperature expansion accounts quite well for the data. 
Evidently the frustration inherent to the kagome AFM prevents much evolution of equal 
time correlations between high and low temperatures in KCr3(OD)6(S04)2. 

V. SUMMARY AND CONCLUSION 

In summary, we have shown that long range AFM order exists below a second order 
transition at TAr=1.8K in the S =3/2 jarosite kagome system KCr3(OD)6(S04)2. The 
ordered structure preserves the chemical unit cell and we have argued that this type of 
structure is favored by interplane coupling. The ordered moment per Cr ion in a sample with 
76% occupancy on the chromium sublattice is 1.1 (3);^^ per occupied chromium site which 
is only one third of that corresponding to Neel order for 5* =3/2 spins. Concomitantly spin 
fluctuations are unusually strong in this system and cannot be accounted for by conventional 
spin wave theory. The inset to Fig. ^ summarizes | < S > \/S versus S for various kagome 
systems studies so far. The data indicate that the weakly diluted S =3/2 kagome AFM 
may be close to a T=0 quantum critical point separating a phase with | < S* > | 7^ from 
a phase with | < 5 > | = . 
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APPENDIX A: NORMALIZING POWDER NEUTRON SCATTERING DATA 

1. Neutron Scattering Cross Sections 

The magnetic neutron scattering cross section from a powder sample at wave vector 



transfer, Q, and energy transfer, huj can be written 

^^TVr^e— «.25(Q..) (Al) 

where rg = 5.38 fm, e"^^'-'^-' is the Debye Waller factor (which we generally take to be unity 
at low Q and low T), and S{Q,uj) is the spherically averaged scattering function: 

d^l^ 1 

The scattering function is given by 

rr' 

where f{Q) is the magnetic form-factor of a single Cr ion [p9[] . 



The Bragg scattering cross section for a powder sample takes the form 

^ = NY.Ar)6{Q-r) (A4) 

where the summation is over reciprocal lattice vectors, r, with different lengths. The Q- 
integrated nuclear scattering intensity is 

X^(r) = £^ E \FN{r')\' (A5) 

r '|=T 
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where v* is the volume of the reciprocal lattice unit cell and the structure factor is 

Fjv(g) = e-^'^'^^b^ e'"^-^ (A6) 

d 

For magnetic Bragg scattering the Q-integrated intensity is 

X^M=ro^£^ E \y^r%\\Fu{r')f-\f.Fu{r')\\ {A7) 

\t '\=t 

where the magnetic vector structure factor is given by 

FniQ) = E ('''^'^'^^ < > e'"^'' (A8) 

d 

2. Normalizing the measured Intensities 

The detector count rate is related to the scattering cross section through an expression 
of the form: 

i{Q, huo) = c{huj) I dQ'hduj'ngUQ -Q\uj~ ^^'), (A9) 

where the spectrometer resolution function is normalized to unity: 

j dQ'hduj'UQ^ {Q-Q',uj-uj') = l (AlO) 

The resolution function and the energy dependence of C{hu)/C{0) can be calculated given 
the spectrometer set up [Q. For experiments with fixed incident energy it is a good ap- 
proximation to use 

where k and k' are the incident and scattered neutron wave vectors, and is the reciprocal 
lattice vector of the Bragg reflection being used as analyzer. The value C(0) depends on 
many factors such as the neutron flux, monochromator and analyzer reflectivities, collimator 
transmissions, and detector efficiencies. We determined C(0) by measuring the Q-integrated 
intensity of nuclear Bragg peaks with known values of X/v(r): 
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With this number it is straight forward to derive absolute values of Tm{t) and S{Q,hLo) 
from the detector count rates associated with magnetic scattering: 

^^mIt) = -^^^j^ J lBr^gg{Q,huj)dQhduj (A13) 
Resolution smearing was neglected in the last expression. 
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TABLES 

TABLE L Positions within space group R3m and occupancies per unit cell of atoms in 
KCr3(OD)6(S04)2 at T=20K. The numbers were derived from the data shown in Fig. 1 by Rietveld 



analysis using GSAS [31|. Isotropic Debye-Waller factors, exp(— < > Q^), were used and the 
resulting overall reduced = 3.8. 
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0.1341(1) 
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0.106(2) 


H 
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0.2218(2) 
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0.077(2) 
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TABLE II. The measured and calculated Q-integrated intensities of Bragg peaks for 
KCr3 (00)6(804)2 powder sample B. Magnetic and nuclear peaks were measured with the same 
experimental configuration. Magnetic Bragg intensities were obtained from the difference between 
the scattering intensity at 0.4K and 3.6K shown in Fig. |^ (a). T^j'(t) was calculated using Eq. 
(A7) for the spin configuration shown in Fig. |2| with {g/2) < Sj- >= 0.55. Comparing measurement 



to calculation we obtain a residual coefficient R = J2r \I°^%t) - Z'"'\t)\/ Y.r'^°^'{T) = 0.17. For 
nuclear Bragg peaks we got R = 0.14 when comparing 

to The latter values were 

calculated from Eq. ([A5|) and Table 1. 



T 

(h,k,l) 


(A-i) 


(10-18 m) 


^S'(r) 
(10-18 ni) 


(10-18 m) 


(10-18 m 


1 1 


1.070 


0.061(6) 


0.0581 


0.32(3) 


0.313 
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1.115 


- 
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2.21(4) 
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1 2 
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0.03(3) 
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1 1 
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0.06(1) 
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10 4 
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2 i, 1 1 3 
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4.683 
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2 7 
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1.710 


2 2 
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FIGURES 

FIG. 1. Powder diffraction data from sample B at T = 20 K taken on the 32 detector powder 
diffractometer BT 1 at NIST. We used tlie Cu (311) monochromator with A = 1.540 A and 
colHmations 15' — 20' — 7'. The Une through the data shows the Rietveld fit which led to the 
structural parameters shown in Table 1. The bottom panel shows the difference between model 
and data. 

FIG. 2. Stacked kagome lattices as they occur in jarosites. A, B, and C denote three mag- 
netic sublattices. (a) is the spin configuration we propose for KCr3(0D)e (804)2 which satisfies 
inter-plane exchange interactions (shaded triangle), (b) Shows that a stacking of kagome layers 
with the V3 X ^3 AFM structure cannot satisfy inter-plane coupling. 

FIG. 3. (a) Temperature dependence of the magnetic Bragg intensity at Q = (012) which is 
proportional to the squared staggered magnetization, (b) Temperature dependence of the specific 
heat divided by temperature for two samples with different Cr concentration in the kagome layers. 
A 100% Cr3+ and B: 75-90% Ci^+. 

FIG. 4. Difference between low and high temperature neutron diffraction data. The data were 
obtained with Ei = 13.7meV, collimations 40'— 36 — 36' and no analyzer, (a) Shows I(0.4K)-I(3.6K) 
probing spin correlations which develop between T=3.6K and T=0.4K. (b) Shows I(0.4K)-I(50K) 
which is a measure of the entire low temperature equal time spin correlation function. The solid 
symbols were obtained by integrating the constant-Q scans shown in Fig. |5[ Dotted and dashed 
lines are described in the text. 
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FIG. 5. Constant-Q scans aX Q = l.lA~^ (filled symbols) and Q = 0.55A~^ (open sym- 
bols). Triangles were obtained with ii^j=5meV and collimations 60 — 40 — 80 — 80 , circles 
with £'i=13.7meV and 60 — 20 — 40 — 40 and rectangles with Ei=14.7meV and collimations 
60 — 20 — 40 — 40 . Backgrounds were subtracted as follows : For huj < 0.7meV we subtracted 
neutron energy gain data. For Hlo > 0.7meV we subtracted data obtained with the analyzer rotated 
10° from its reflection condition. The inset shows \ < S > \/ S versus S for various kagome systems. 
Squares: S = 5/2 Fe-jarosite systems with long range order [111]; Open circles: frozen moment in 



the spin glass phase of SrCr9pGai2-9pOi9 with p=0.92 and S = 3/2 [15|; Filled circle: long range 
ordered moment in the S = 3/2 Cr-jarosite KCr3(OD)6(S04)2 with 76% occupancy of the kagome 
lattice. 

FIG. 6. Constant-?iw scans at 0.5, 1.1, and 2 meV. The ?itiJ=0.5meV scan was obtained with 
Ei=5meV and collimations 60' — 40' — 80' — 80'. The huj=l.lmeV and huj=2.0meV scans were 
obtained with £'j=13.7meV and collimations 60' — 20' — 40' — 40'. The background subtraction 
scheme is described in the caption to Fig. ^. Solid lines are a single parameter fit to a spin wave 
theory described as in the text. 
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